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Abstract. The i?-unistochastic matrices are a special class of symmetric bis- 
tochastic matrices obtained by taking the square of the absolute value of each 
entry of a Hermitian unitary matrix. We examine the geometric relationship 
of the convex hull of the n by n if-unistochastic matrices relative to the larger 
convex set of n by n symmetric bistochastic matrices. We show that any line 
segment in the convex set of the n by n symmetric bistochastic matrices which 
passes through the centroid of this convex set must spend at least two-thirds 
of its length in the convex hull of the n by n .ff-unistochastic matrices when 
n is either three or four and we prove a partial result for higher n. A class of 
completely positive linear maps called the self-dual doubly stochastic maps is 
useful for studying this problem. Some results on self-dual doubly stochastic 
maps arc given including a self-dual version of the Laudau-Streater theorem. 



1. Some Convex Subsets of the Birkhoff Polytope 

A doubly stochastic (or bistochastic) matrix is a square non-negative matrix each 
of whose rows and whose columns sum to one. We let B n denote the set of n by n 
doubly stochastic matrices, B n is sometimes called the Birkhoff polytope. We will 
consider certain subsets of doubly stochastic matrices. The set of n by n symmetric 
bistochastic matrices, denoted B^ m is the set of all symmetric matrices in B n . If 
U is a n by n unitary matrix, let Uo U be the n by n matrix whose (i, j)th entry is 
\uij\ . U o U is clearly doubly stochastic; any doubly stochastic matrix that can be 
obtained in this way is called unistochastic. (We follow standard notation is letting 
A o B denote the Schur product of two n by n matrices A and B; recall that A o B 
is the n by n matrix whose (i,j)th entry is aijbij). If U is a Hermitian unitary 
matrix, then U o U is a symmetric bistochastic matrix. Any matrix of the form 
U o U where U is a Hermitian unitary matrix is called an ii-unistochastic matrix, 
if-unistochastic matrices were defined and first studied in |17j . When n > 3, 
there exist doubly stochastic matrices which are not unistochastic and symmetric 
bistochastic matrices which are not .ff-unistochastic. 

/° i i\ 

Example 1. The matrix ^ o \ \ is an example of a matrix which is symmetric 

bistochastic but not unistochastic. 

B n and B^> m are both convex sets. It is a well known result of Birkhoff that 
the extreme points of B n are the permutation matrices. Katz [7] proved that all 
extreme points of B s ^ m are of the form i(P + P T ) where P is a permutation matrix 



2010 Mathematics Subject Classification. 15B51, 15B48. 

Key words and phrases, doubly stochastic matrices, H-unistochastic matrices, completely pos- 
itive linear maps, correlation matrices. 

1 



2 



COREY O'MEARA AND RAJESH PEREIRA 



whose corresponding permutation contains no even cycle of length greater than or 
equal to four. We state Katz's theorem more exactly below. 

Theorem 1. [7] A matrix M is an extreme point of the convex set of n by n 
symmetric bistochastic matrices if and only if there exists a permutation matrix P 
such that PMP T is a block diagonal matrix with diagonal blocks {-Bn fc }fcLi w ^ ere 
to > 1, nk is either odd or equal to two for each k, Y^k=i nk ~ 71 an d w here B nk is 
therik byrik matrices defined as follows: B\ = (1), T$2 = (? o) and for nk > 3, B nk 
is the matrix whose (i,j) entry is \ if i ~ 3 = il ( mod nk) and is zero otherwise. 

We note that B 3 is the matrix given in example [1] In fact for any odd n greater 
than or equal to three, B n is a symmetric bistochastic matrix which is not unis- 
tochastic. We will call the decomposition of an extreme point of B s ^ m given in the 
above theorem the Katz decomposition. 

Since every permutation matrix is clearly unistochastic, the set of n by n doubly 
stochastic matrices is the convex hull of the set of n by n unistochastic matrices. 
When n > 3, the convex hull of the n by n the convex hull of the iJ-unistochastic 
matrices is strictly smaller than the set of symmetric bistochastic matrices, since 
any extreme point of B^ m whose Katz decomposition has blocks of size three by 
three or larger will not be unistochastic. 

We will let B^ denote the convex hull of the n by n i/-unistochastic matrices. We 
note that all three convex sets B n , B^ m and B\ are all invariant under permutation 
similarities (i.e. maps of the form X — > PXP T ). The centroid of each of these 
convex sets is clearly the n by n matrix all of whose entries are ^ as it is the unique 
point in B n which is fixed by all permutation similarities. We will follow |17j and 
let W n denote the n by n matrix all of whose entries are —. We will examine 
the relationship between the convex sets B^ m and B\. We then apply the theory 
of completely positive linear maps to some of our results; we review the basics of 
completely positive linear maps in the next section. 

2. Introduction to Completely Positive Linear Maps 
In this section we review the basic theory of completely positive linear maps. 

Definition 1. [TB] A linear map $ : M n (C) M n (C) is 

(1) positive if (&(A) is positive semidefinite whenever A is positive semidefinite. 

(2) completely positive if $ <S> Ik is a positive linear map from A / Ifc„(C) — > 
M kn (<C) for allkeN. 

(3) trace-preserving if Tr ((&(A)) = Tr(A), for all A € A4 n . 

(4) unital if <$>{I n ) = I n . 

We will call a unital completely positive trace preserving linear map a doubly 
stochastic map. 

The following theorem of Choi provides a necessary and sufficient condition on 
the decomposition of linear map to be completely positive. 

Theorem 2. [3] v4 linear map $ : A4 n (<C) — >■ M n (C) is completely positive if and 
only if it admits the following form 

r 

(1) ${A) =Y,ViAV* 

i=l 

where {V,} G M n (C). 
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The matrices {V{\ are called the Kraus operators for a given map <£>, and the 
choice of such Kraus operators of $ is non-unique. The number of terms in the 
summation will be at a minimum whenever the Kraus operators are chosen to be 
linearly independent, this minimal value of r < n 2 is called the rank of $. The 
rank one doubly stochastic maps arc precisely the unitary transforms. 

Definition 2. Let <I> : A4 n (C) — > Ai n (C) be a completely positive map. & is a 
mixed unitary map if there exists a collection of unitary matrices {Ui} i=1 such that 

in 

$(A) =J2PiUiAU* 

i=l 

where Y^lLiPi — 1 an d Pi > 0- 

Mixed unitary maps are also often called random unitary maps in the literature; 
in using the term mixed unitary we are following [19] . Landau and Streater [10] 
have shown that every doubly stochastic map on M.2{ < C) is a mixed unitary map 
but there exist doubly stochastic maps on .A/f„(C) for all n > 3 which are not mixed 
unitary. Much more about doubly stochastic and mixed unitary maps and their 
application to physics can be found in [15] . 

3. Geometrical Relations between some convex subsets of the 

BlRKHOFF POLYTOPE 

The following general problem has been considered in many different contexts. 

Problem 1. Let V be a finite dimensional real of complex vector space and let A 
and B be two compact convex subsets of V with B C A. Suppose A and B have 
a common centroid c G B . Find X, the largest constant in [0, 1] for which the set 
inclusion XA + (1 — A)c C B is satisfied. 

Some norm inequalities are in fact essentially the solution to this problem where 
A and B are the unit balls of the norms in question. This problem has also been 
explored in the context of doubly stochastic maps and correlation matrices. In [19] . 
the case where A is the set of doubly stochastic maps on M„(C) and B is the set 
of mixed unitary maps on M n (C) was considered. Here a lower bound of n 2 1 _ 1 was 
found for A for each fixed n. In |14j . the case where A is the n by n real correlation 
matrices and where B is the the convex hull of the n by n real rank one correlation 
matrices was considered. It was shown that A = | when 3 < n < 6 and bounds were 
given for higher dimensional cases. The motivations for the latter two problems are 
quantum information and semidefinite programming respectively. 

We will consider the case of problem [T] where A is the set of n by n symmet- 
ric bistochastic matrices and B is the convex hull of the n by n If-unistochastic 
matrices. 

Theorem 3. Let X n be the largest constant in [0, 1] for which X n M + (1 — X n )W n £E 
for all n by n symmetric bistochastic matrices M. Then Ai = A2 = 1, A3 = 
A 4 = §. 

That Ai = A2 = 1 is an immediate consequence of the fact that all 1 by 1 and 2 
by 2 symmetric bistochastic matrices are iJ-unistochastic. We split the remainder 
of this result into two parts as the proofs are quite different for each case: 

Proposition 1. A3 = |. 
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Proof: We note that there is only one extreme point of B^ ym which is not H- 

/o | i\ 

unistochastic; this is the matrix M = \ o \ from example [T] So we only need 
show that kM + (1 — k)W% is in ff 3 for fc = | but not in H 3 if fc > §. Note 

(i * * \ / _i 2 3. \ 

9 9 9 \ / 3 3 3 \ 

llf is £f-unistochastic since I - 1 I is a Hcrmitian 
999/ \ f f ~ 3 / 

unitary matrix. 

We note that if n is odd and U is an n by n Hcrmitian unitary matrix, then 
|fr([/)| > 1 since all the eigenvalues of U are plus or minus one. Hence by the 
Cauchy-Schwarz inequality, X)"=i \ u a\ 2 — n \ S"=i u a\ 2 ^ ^ which implies that the 
trace of any matrix in B% must be greater than or equal to ^ when n is odd. Since 
tr(kM + (1 - fc)W 3 ) = 1 - fc < 5 when fc > | our result follows. ■ 



Proposition 2. A4 = |. 



Proof: We note that up to permutation similarity there is only one extreme point 

'0 i i N 



of Bf 1 "" which is not H- unistochastic; this is the matrix M = 2 2 I . Now 

2 5 1 

,6 6 01, 



' 







let X = 

1 1 n 1 I • This matrix is arithmetic mean of the three Hcrmitian 



333" 



permutation matrices corresponding to the three permutations (12)(34), (13)(24) 

(ill n\ 
9 9 9 u \ 
Iff I ; it follows from the proof of 
9 9 9 °/ 
0001/ 

the previous proposition that Y e B\. We note that §M + \W± = jX + jY 6 B\. 
Hence A4 > |. 

We now prove that if A e B\ , then 3an + 3a 2 2 + 3a33 — 044 > 0. We only need 
show this for the special case where A is ff-unistochastic; the general case where A 
is in the convex hull of ii-unistochastic matrices follows from convexity. Let U be a 
Hcrmitian unitary matrix such that U o U = A. Note that tr(U) is an even integer 
and U44 € [-1, 1], hence \uu\ < |"44 - tr(U)\ = \un + u 2 2 + u 33 \ < J2i=i \ u u\- 
Squaring both sides of the last inequality and then using Cauchy-Schwarz, we get 
a 44 = |u 44 | 2 < (X) i= i Wu\) 2 < 3X)i=i \ u ii\ 2 = 3 X)Li a "- Now tnc diagonal entries 
of kM + (1 - k)Wi arc i=^, i=^, ^). Hence if fcM + (1 - fc)W 4 G Bj, then 

9 — 9fc > 1 + 3fc which means that fc < |. Therefore A 4 = |. ■ 

We note that the number | appears both as a solution to the small dimensional 
cases of both the correlation matrix problem from [T3] and our problem on sym- 
metric bistochastic matrices. While this may be coincidental, we will soon derive a 
result which suggests some connections between these two problems. We first need 
to introduce the class of self-dual completely positive linear maps and explore its 
connections to both the symmetric bistochastic matrices and the real correlation 
matrices. Wc do this in the next section, which culminates in a partial result of our 
problem for general n. The exact value of X n for n > 5 remains an open problem. 
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4. Self-dual linear maps and correlation matrices 



As the names suggest, there are very close connections between doubly stochastic 
matrices and doubly stochastic maps. (Some of these connections will be reviewed 
later on in the paper). The doubly stochastic maps are in a sense an analogue of 
the doubly stochastic matrices. We would like to find a similar analogue of the 
symmetric bistochastic matrices. We note that Ai n (C) is a Hilbert space under the 
inner product (A, B) = Tr(AB*). The adjoint or dual to the CP map $ is the linear 
map $* which satisfies (A, $(-B)) = B) for all A, B. It follows from simple 

calculation that the dual of = £T =1 V t AV* is = ££ =1 V* AV t . Hence 

the dual of a CP map is a CP map. We say that a CP map is self-dual if $ = $*. 
The self-dual doubly stochastic maps are the analogue of the symmetric bistochastic 
matrices. We study the properties of self-dual completely positive linear maps in 
this section. 

The following theorem is straightforward and can be found in j!Zu] . 

Theorem 4. [20] A completely positive map $ : A4 n (C) — > A4 n {C) is self-dual if 
and only if there exists a set of n by n Hermitian matrices {Hi} such that: 



for all A e M n (C). 

Proof: The if direction is immediate. Now suppose $>(A) = Y^i=i ViAV* is self- 
dual. Since <£> = $*, the self-dual map $ can be expressed as 4> = . Since any 
matrix V can be decomposed via Cartesian decomposition as V = K + iL, where 
K and L are n x n Hermitian matrices, decomposing each Kraus operator in this 
manner yields the result &(A) = Y^l=i HiAKi + LiAL t ■ 

We recall the result that two sets of Kraus operators {Vi} and {Wj} represent 
the same completely positive map $ if and only if there exists a unitary matrix 
U = (uij) such that V = J^j u ijWj [!]■ While not every set of Kraus operators 
of a self-dual map will be Hermitian, theorem [4] guarantees that there is always at 
least one set of Kraus operators which are Hermitian for a self-dual map. We now 
introduce a self-dual version of the concept of a mixed unitary map. 

Definition 3. Let $ : M. n (C) — > A^ n (C) be a self-dual completely positive map. 
$ is a mixed Hermitian unitary map if there exists a collection of Hermitian unitary 
matrices {Ui}™ =1 such that 



where Y^iLiPi = 1 and pi > 0. 

The relationship between the mixed Hermitian unitary maps and the self-dual 
doubly stochastic maps can be thought of as being analogous to the relationship 
between the mixed unitary maps and the doubly stochastic maps. There is one key 
difference which we will describe in the next paragraph. 

Let $ be a doubly stochastic map on A4 n (C), and define A$(D) = /„ o 
for all n by n diagonal matrices D. Since the n by n diagonal matrices form an 
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n dimensional vector space with the obvious standard basis, we will identify A$ 
with its n by n matrix representation. Hence if $ is a doubly stochastic map, then 
A$ is a doubly stochastic matrix. If <3> is a self-dual doubly stochastic map, then 
A$ is a symmetric bistochastic matrix. Conversely, if M is a doubly stochastic 
matrix (resp. symmetric bistochastic matrix), there exists a doubly stochastic map 
(resp. self-dual doubly stochastic map) $ such that A<[, = M. The converse 
statements follows from considering the rank one map = PAP T where P is 

any permutation matrix. Then it can be shown that A$ = P. By linearity, maps of 
rank r < n 2 must be doubly stochastic. If U is a unitary matrix and $>(A) = U AU*, 
then A$ = U o U; hence by linearity of A we have the following result: 

Proposition 3. Let $ be a completely positive linear map on M n (C). //$ is mixed 
Hermitian unitary, then A$ lies in the convex hull of the n by n H -unistochastic 
matrices. 

This gives us a key difference between the self-dual case and the general case. 
Given a doubly stochastic linear map $, the matrix A$ will never give us any 
information as to whether $ is mixed unitary since the convex hull of the n by n 
unistochastic matrices is the entire set of n by n doubly stochastic matrices. 

We introduce the set of correlation matrices which will be useful in constructing 
doubly stochastic maps. 

Definition 4. Let C £ A4 ra (C) be a positive semidefinite matrix with cu — 1 for 
all 1 < i < n. Then we say that C is a correlation matrix. 

For any n by n correlation matrix C, we follow the convention in |12| by letting 
$c denote the completely positive map which maps A £ M n to the Schur product 
A o C. We note that if C is an n by n correlation matrix, then $c is a doubly 
stochastic map with $ c t being its dual. Hence &(A) = A o C is self-dual if and 
only if C is a real correlation matrix. Any completely positive linear map which 
leaves every diagonal matrix unchanged must be of this type; this statement of the 
result was given independently in [5] and [T2J. See for instance proposition 2.1 of 
Kye [H] and the discussion previous to it or theorem 1 from [T2]. A more general 
result was found even earlier in |T5] and was rediscovered in [2]. We will not use 
the more general result here and refer to [21 El HI] for details. 

Proposition 4. [SJ 112) Let $ be a completely positive linear map on the n by n 

matrices such that = D for all diagonal matrices D . Then there exists an n 

by n correlation matrix C such that $(^4) = $c{A) — C o A for all A € M n (C). 

The next two results provide conditions on whether such maps are mixed unitary 
or mixed Hermitian unitary. 

Theorem 5. Let C be a correlation matrix. Then $c is a mixed unitary map if 
and only if the associated correlation matrix C is a convex combination of rank 1 
correlation matrices. 

Proof: Suppose &(A) is a mixed unitary map, then $(A) — Y^k=iPkUkAU^ where 
{Uk} r k=1 are unitary, X)fc=iPfe = 1 ana - Pk > 0. Let Ei be the matrix with a one in 
its ith row, ith column entry and zeros everywhere else. Since $(i?;) = Ei which 
is rank one, UkEiU^ = Ei for 1 < i < n and 1 < k < n. Therefore Uk commutes 
with all diagonal matrices which means it must be a diagonal matrix. Let be the 
vector consisting of the diagonal elements of Uk, then VkV% is a rank-one correlation 
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matrix and C = J2k=i PkVkV^. Similarly if C = J2k=i PkVkV^ where VkV^. are rank- 
one correlation matrices, then we let Uk be the diagonal unitary matrices formed 
using the entries of Vk and hence C o A ~ Y^k=i PkUkAU^. ■ 
By using the same argument as above, we obtain: 

Theorem 6. Let C be a real correlation matrix. Then $c * s o. mixed Hermitian 
unitary map if and only if the associated correlation matrix C is a convex combi- 
nation of real rank 1 correlation matrices. 

In light of this, it would be useful to be able to test whether a real correlation 
matrix C is the convex hull of the real rank one correlation matrices. The following 
result is a very slight modification of a result from [1]. 

Theorem 7. Let C be an n by n real correlation matrix and let P n be the collection 
of all subsets of {1,2, ...,n}. Then C is in the convex hull of the real rank one 
correlation matrices if and only if there exists a function f : P n — > R with the 
following properties: 

(1) /(0) = i 

(2) f({i,j}) = Cij forl<i< j<n 

(3) £T e p„M) |SnT| /Cn > for all S E P n 

Proof: Let K be the set of all functions / : P n — > R which satisfy conditions one 
and three from the statement of the theorem above. It is clear that K is a convex 
set. For all A g P n , let f A (T) = (-l)l AnT l. Since £tgp„ (-l) |SnT| /A(T) = 
except when S = A, {/a}a£P„ are exactly the extreme points of K. Let Ca be the 
n by n matrix whose main diagonal elements are all one and whose (i,j)th entry 
is j}) for all i ^ j. Since {CUIagP™ is the set of all n by n real rank one 

correlation matrices, the result follows. ■ 

We can use Theorem [7] to give us a lower bound for a rank-dependent variant of 
the problem from |14j . 

Corollary 8. Let C be an n by n real correlation matrix and let r = rank{C). 
Then -C + £=1/ is in the convex hull of real rank one correlation matrices. 

Proof: There exists a real r by n matrix A such that C = A T A. We note that all 
entries of A must have absolute value less than or equal to one. Now we construct 
the following function g : P n — > R. Let g(0) = r and let g(T) = Y^l=i YijeT a «i 
whenever T is non-empty. We note that g{{i,j}) = c^ and for any S £ P n . 

ErepJ-^ SnTl 9(T) = £U EtcpJ-I) 1 ^ 1 U jeT *U = EU II"=i(l+(-l) l{j ' } 
0. Now take the function f(T) = r 1 g(T) and apply / to theorem [7] to show that 
-C + ^-^r-I is in the convex hull of real rank one correlation matrices. 

The set of all n by n correlation matrices in M. n {¥) over a field F which is either 
the real or the complex numbers form a compact convex set that is usually denoted 
e„(F). As with most convex sets, there has been interest in the properties of its 
extreme points. It was shown by Christensen and Vesterstrom [5] and Loewy |13j 
that e n (C) has extreme points of rank r if and only if r < y/n. Later, Grone, Pierce 
and Watkins [B] found an analogous result for the set e„(R) which they showed 
has extreme points of rank r if and only if r(r + 1) < 2n. A simple derivation of 
these results along with an efficient algorithm for checking whether a correlation 
matrix is extremal was given by Li and Tarn in The Li- Tarn characterization 
of extremal correlation matrices has been generalized to the infinite-dimensional 
case by Kiukas and Pellonpaa in [8] . 



8 



COREY O'MEARA AND RAJESH PEREIRA 



We let p n ^ r be the smallest number positive number less or equal to one than one 
for which p n , r C + (1 — p n ,r)In is in the convex hull of rank one correlation matrices 
whenever C is an n by n real correlation matrix of rank at most r. From |14j . we 
know that ,03.2 = P4,2 = P5.2 = P&,2 = P6,3 = §■ We also know that when n > 3 
and fc>2,i<p nr <| where the lower bound follows from corollary [5] and the 
upper bound follows from the fact that {p n ,r} is a nonincreasing sequence in both 
n and r. 

We now return to our main problem. First we recall some notation. For any 
natural number n, let F n denote the n by n Fourier matrix which is the unitary 
matrix whose [j, fcjth entry is n » l . If U is any n by n unitary, then 

denotes the unitary conjugation map which sends A £ M n to UAU* . 

Proposition 5. Let n > 3, let q <E N with mq = n and let M be a n by n be an 

extremal symmetric bistochastic matrix having q Katz blocks all of which are m by 
to. Then p n , 2 M + (1 - p„, 2 )W„ £ B h n . 

Proof: In this proof we note that we are using the symbol o to denote composition 
and not Schur multiplication. Let n = mq and let M be the extremal symmetric 

27ri 

bistochastic matrix having q Katz blocks all of which are m by m. Let u> = e~ 
and q = ^. Then M is permutationally similar to B m <£> I q , so without loss of 
generality let us assume M = B m <E> I q . 

Let C m be the to by to correlation matrix whose (J, fc)th entry is cos( 27r ^ fc ) ) = 
— — y 1 — -. We note that the (j, fc)th entry of C m is equal to (e 3 -,efc) where eu = 
(cos(^ 2 ^), sin(^^)) T is the unit vector with polar angle radians for all k; hence 
C m is a rank two correlation matrix. Let v p be the pth column of F m , then 
is an to by m matrix whose (j, fc)th entry is ic;(P _:L )w _fc ). Hence ^C ra (%>Wp) = 
5 (i>p— ii>p_i + Up+ifp.).!) where the additions p — 1 and p + 1 are taken mod p. Let 
-Ep denote the m by m matrix with a one in its (p,p)th entry and zeros everywhere 
else. Note that v p v* = F m E p F m . Let $ = Tp^ o $ Cm r Fp , then § n (E p ) = 
7}(E p -i + E p+ i) where the additions p — 1 and p + 1 are taken mod to. 

Note that A^, m(Stdq = M. Let U = F m <E)F q , then <S> m ®id q = Tu*°®qC m ®w q Tu. 
Let S„ = Tjj- o $j n o ry. We note that for any n by n diagonal matrix D, S(£)) = 
M- ) j n anc [ so = W n . Since <7C m (8 is an n by n rank two correlation 
matrix, B = p n ^{<iC m ® W q ) + (1 — p n ,2)In is a convex combination of rank one 
matrices. Hence by theorem[6l \& = p n , lif&m ® id q ) + (1 — Pnfi^n = IV» o $ B oTy 
is mixed Hermitian unitary map. Therefore p n , 2 M + (1 — p n! 2)W n = A* € Bl[ by 
proposition [3] 

5. The self-dual Landau-Streater theorem 

The set of doubly stochastic maps on A4 n (C) is convex and as with most con- 
vex sets, there has been interest in finding the extreme points. Since the unitary 
transforms Tu(A) = UAU* are the rank one doubly stochastic maps, they must be 
extremal. Are they the only extreme points? Landau and Streater |10j have shown 
that every extreme doubly stochastic map on A^C) must be a unitary transform, 
but there exist extremal doubly stochastic maps on M n (C) which are not unitary 
transforms for all n greater than or equal to three. It follows from this that every 
doubly stochastic map on M. 2 { < £) is a mixed unitary map but there exist doubly 
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stochastic maps on Af n (C) for all n > 3 which are not mixed unitary. In this sec- 
tion, we show that the situation for self-dual maps is similar; every self-dual doubly 
stochastic map on A^2(C) is a mixed Hcrmitian unitary map but there exists doubly 
stochastic maps on Af„(C) for all n > 3 which are not mixed Hcrmitian unitary. 
Our results on correlation matrices give us a particularly simple construction of the 
latter. 

( 1 * *\ 

Example 2. Let C = I -\ i — § . C is a real correlation matrix which is in 

the convex hull of complex rank-one correlation matrices but not in the convex hull 
of the real rank-one correlation matrices. (In fact, C is a rank 2 extreme "point of 
£3(R) / ). Hence the map A — > AoC is a self-dual mixed unitary map which is not 
mixed Hcrmitian unitary. 

Since the convex set of n by n real correlation matrices has rank two extreme 
points for all n > 3, there exist self-dual doubly stochastic maps on M. n (C) which 
are not mixed Hcrmitian unitary. We can also use the fact that there are symmetric 
bistochastic matrices which are not the convex combination of iJ-unistochastic 
matrices to give another class of examples. 

Example 3. Let n > 3 and let P be an n by n permatation matrix corresponing 
to the permutation (123). Let = \(PAP T + P T AP). Then consists of 

the extremal symmetric bistochastic matrix which is the direct sum of a single Katz 
block i?3 with n — 3 Katz blocks B\. Since A$ is not in the convex hull of the H- 
unistochastic matrices, <j> cannot be mixed Hermitian unitary. Note however that 
this map like that of example^ is both self-dual and mixed unitary. 

We now prove that for n = 2, all the self-dual doubly stochastic maps are mixed 
Hermitian unitary. This result is the self-dual analogue of the Landau-Streater 
theorem. 

Theorem 9. Let <1> : A^C) — > A^C) be a doubly stochastic map. Then <f> is 
self-dual if and only if it admits a Kraus decomposition of the form 

m 

${A) =^p i R i AR i 

i 

where {Ri}™=i is a collection of unitary Hermitian matrices and X)i=iP« = 1 
Pt > 0. 

Proof: The if direction is immediate. We prove the only if part. Let ^ be an 
extreme doubly stochastic map on A^CC). Then wc will show that the map $ = 
£^ admits a Kraus decomposition composed of Hermitian unitary matrices and 
the general case will follow by linearity. Since the extreme points are exactly the 
unitary transforms, let U = VDV* be the spectral decomposition of the single 
Kraus operator U such that = ^o* fl o 4>*(A), where 4> '■ A — > VAV* 

and ^ d '■ A — > dad +d ad ^ Clearly, the diagonal matrix resulting from the 
spectral decomposition is D = diag^e* 9 , e _i ^). By Proposition [4j is a Schur 
map expressed as 

T / a \ 1 cos (8 + 6) 1 . 
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Since every 2x2 real correlation matrix is a convex combination of real rank 1 
matrices, it is clear that = cj> o * o <j>*{A) = J2LiP*( v ~D t V*)A(V D t V*) = 

X)|=i PiRiARi where diag(D,) = vi such that Wi = VivJ '. ■ 

6. Conclusion & Outlook 

We conclude this paper by summarizing our main problem and its possible con- 
nection with the correlation matrix problem studied in |14j . We summarize the 
relationships between the mathematical objects studied in the paper and formulate 
a conjecture based on our partial results. 

While the convex hull of the n by n unistochastic matrices is the set of n by 
n doubly stochastic matrices, B% the convex hull of the n by n if-unistochastic 
matrices is a proper subset of B n ym the set of symmetric bistochastic matrices. We 
considered the problem of finding the smallest nonnegativc number A n for which 
the set \ n B%> m + (1 - X n )W n is contained in the set B s n ym . Clearly Ai = A 2 = 1. 
We found the first two non- trivial cases of A„ by showing that A3 = A4 = | . 

In [14], a similar problem was considered. Let e n (M.) be the convex set of all n by 
n real correlation matrices and let Q n be the convex hull of the set of n by n real rank 
one correlation matrices. The problem in [14) is to find the smallest nonnegativc 
number p n for which the set p n e n (BL) + (1 — p n )I n is completely contained in Q n . 
The known values of p n are pi = pi = 1 and P3 = Pa = P5 = p& = ^ [H] • We note 
that A„ = p n when n < 4 which accounts for all of the known values of A„ . There 
are also other connections between these two problems. Specifically, $c(^4) = Co A 
is a self-dual doubly stochastic map if and only if C is a real correlation matrix 
and $c is mixed Hermitian unitary if and only if C is in Q n . Furthermore, if $ is 
any self-dual doubly stochastic map, then A$ is a symmetric stochastic matrix; if 
in addition $ is mixed Hermitian unitary then A$ € £>'>. These two results were 
used to prove proposition [5] which also suggests a connection between A„ and p n . 
Motivated by this, we offer the following conjecture: 

Conjecture 1. A„ = p n for all n E N. 

In the previous section, as a further application of these concepts, we show that 
every self-dual doubly stochastic maps on M n (C) is mixed Hermitian unitary if and 
only if n < 2. This is the self-dual version of the well-known result of Landau and 
Streater that every doubly stochastic maps on M n (C) is mixed unitary if and only 
if n < 2. 
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